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Abstract. This talk summarizes the present status of a program to quantitatively relate data 
from the Relativistic Heavy Ion Collider (RHIC) on collective expansion flow to the Equation 
of State (EOS) of hot and dense strongly interacting matter, including the quark-gluon plasma 
and the quark-hadron phase transition. The limits reached with the present state of the art and 
the next steps required to make further progress will both be discussed. 



1. The QCD Equation of State and ideal fluid dynamics 

With relativistic heavy-ion collisions one explores the phase diagram of strongly interacting bulk 
matter in the regime of high energy density and temperature. Lattice QCD (LQCD) tells us 
that for zero net baryon density QCD matter undergoes a phase transition at T cr = 173 ±15 MeV 
from a color-confined hadron resonance gas (HG) to a color-deconfined quark-gluon plasma 
(QGP). The critical energy density e cr ~ 0.7 GeV/fm 3 corresponds roughly to that in the 
center of a proton. At the phase transition, the normalized energy density e/T rises rapidly by 
about an order of magnitude over a narrow temperature interval AT < 15 — 20 MeV, whereas 
the pressure p/T 4 (which is proportional to the grand canonical thermodynamic potential) is 
continuous and rises more gradually (Fig. [J). Both seem to saturate at about 80-85% of the 
Stefan-Boltzmann value for an ideal gas of noninteracting quarks and gluons, the energy density 
more quickly (at about 1.2 T cr ), the pressure more slowly. Above about 2T cr , the lattice data 
follow the Equation of State of an ideal gas of massless particles, e = 3p. 

For many years this observation has been interpreted as lattice QCD support for the 
hypothesis of a weakly interacting, perturbative QGP. The recent RHIC data taught us (as 
I will show) that this interpretation was quite wrong. 

It was recognized over 3 decades ago (see review [2]) that information about the EOS of 
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Figure 1. The normalized energy density e/T 4 (left) and pressure p/T 4 (right) from lattice QCD pQ 
for 0, 2 and 3 light quark flavors, as well as for 2 light + 1 heavier (strange) quark flavors. Horizontal 
arrows on the right indicate the corresponding Stefan-Boltzmann values for a non-interacting quark-gluon 
gas. 
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strongly interacting matter can be extracted by studying the collective dynamics of relativistic 
heavy-ion collisions. This connection is particularly direct in the framework of ideal fluid 
dynamics which becomes applicable if the matter formed in the collision approaches local 
thermal equilibrium. The latter requires sufficiently strong interactions in the medium that 
local relaxation time scales are shorter than the macroscopic evolution time scale. In this limit 
the local conservation laws for the baryon number, energy and momentum currents, d^jg(x) = 
and dfjT^ u = 0, can be rewritten as the relativistc Euler equations for ideal fluid motion: 

riB = -ns(d-u), e = — (e + p) (d ■ u), 

e + p 1 + cj \eoJ 

The dot denotes the time derivative in the local fluid rest frame (/ = u- df) and V^ 1 the gradient 
in the directions tranverse to the fluid 4-velocity u M . The first line describes the dilution of 
baryon and energy density due to the local expansion rate d ■ u, which itself is driven according 
to (|2j) by the pressure or energy density gradients providing the fluid acceleration. The absolute 
value of the energy density e is locally irrelevant: the initial maximal energy density eo only 
matters by setting the overall time scale between the beginning of hydrodynamic expansion and 
final decoupling, thereby controlling how much flow can develop globally. The details of the flow 
pattern are thus entirely controlled by the temperature dependent speed of sound c 2 s = 

According to the LQCD data, the latter is c 2 s ~ | for T > 2 T cr , then drops steeply near T ~ T cr 
to values near c 2 k, ^ (the "softest point", see Fig. 11 in IdJ), before rising again in the hadron 
resonance gas phase to c 2 « 0.15 [Hj. A key goal of flow studies in relativistic heavy ion collisions 
is to find traces of this "softest point" in the data. 

2. "Flavors" of transverse flow in heavy ion collisions 

Experimentally one studies flow by analyzing the transverse momentum spectra of the emitted 
hadrons. In central (6 = 0) collisions between spherical nuclei, the flow is azimuthally symmetric 
about the beam axis. This "radial flow" integrates over the entire pressure history of the 
collision from initial thermalization to final decoupling ( "freeze-out" ) , due to persistent pressure 
gradients. In noncentral (6^0) collisions, or central collisions between deformed nuclei such as 
uranium the nuclear reaction zone is spatially deformed, and anisotropies of the transverse 
pressure gradients result in transverse flow anisotropies. These can be quantified by Fourier 
expanding the measured final momentum spectrum dN/(dy pj_dpj_d(j)p) with respect to the 
azimuthal angle <p p . For collisions between equal nuclei, the first non- vanishing Fourier coefficient 
at midrapidity is the elliptic flow V2 (p± , b) . Since V2 is driven by pressure anisotropies and the 
spatial deformation of the reaction zone creating such anisotropies quickly decreases as time 
proceeds, the elliptic flow is sensitive to the EOS only during the early expansion stage |5| (the 
first ~ 5 fm/c in semicentral Au+Au collisions [S]), until the spatial deformation has disappeared. 

Depending on the initial energy density (i.e. beam energy), the hot expanding fireball spends 
this crucial time either entirely in the QGP phase, or mostly near the quark-hadron phase 
transition, or predominantly in the hadron resonance gas phase [H], thereby probing different 
effective values of the sound speed c s . To the extent that ideal fluid dynamics is valid in all these 
cases, an excitation function of the elliptic flow V2 should thus allow to map the temperature 
dependence of the speed of sound and identify the quark-hadron phase transition, via a minimum 
in the function V2(y/s) This will be further discussed below (see Section f5, 21 and Fig. EJ). 

3. Model parameters and predictive power of hydrodynamics 

The hydrodynamic model requires initial conditions at the earliest time at which the assumption 
of local thermal equilibrium is applicable, and a "freeze-out prescription" at the end when 
the system becomes too dilute to maintain local thermal equilibrium. Both are described in 
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detail elsewhere jjj. Different approaches to freeze-out invoke either the Cooper-Frye algorithm 
(Hj (used by us), in which chemical freeze-out of the hadron abundances at T cr [2j must be 
implemented by hand by introducing non-equilibrium chemical potentials below T cr |10[ I11 | IT2]. 
or a hybrid approach |13[ 11 4 j that switches from a hydrodynamic description to a microscopic 
hadron cascade at the quark-hadron transition, letting the cascade handle the chemical and 
thermal freeze-out kinetics. While the radial flow patterns from the two freeze-out procedures 
don't differ much, the elliptic flow can be quite different if the spatial deformation of the fireball 
is still significant during the hadronic stage of the expansion, as I will discuss in Sec. 15.21 

We have solved the relativistic equations for ideal hydrodynamics under the simplifying 
assumption of boost-invariant longitudinal expansion (see for details). This is adequate 

near midrapidity (the region which most RHIC experiments cover best), but not sufficient to 
describe the rapidity distribution of emitted hadrons and of their transverse flow pattern which 
require a (3+l)-dimensional hydrodynamic code such as the one by Hirano |lUj . 

The initial and final conditions for the hydrodynamic evolution are fixed by fitting the 
pion and proton spectra at midrapidity in central (6 = 0) collisions; additionally, we use the 
centrality dependence of the total charged multiplicity dN^/dy. I stress that this is the 
only information used from 6^0 collisions, and it is necessary to fix the ratio of soft to 
hard collision processes in the initial entropy production. An upper limit for the initial ther- 
malization time To<0.6fm/c, the initial entropy density so = 110fm~ 3 in the fireball center 
(corresponding to an initial maximal energy density eo~30GeV/fm 3 and an initial central 
fireball temperature To«360MeV«2T cr ), the baryon to entropy ratio, and the freeze-out 
energy density edec = 0.075 GeV/fm 3 are all fixed from the 6 = pion and proton spectra (these 
numbers refer to 200^4GeV Au+Au collisions at RHIC |12j). The initial entropy density profile 
is calculated for all 6 from the collision geometry, using the Glauber model with soft/hard ratio 
as fixed above. Except for pions and protons, all other hadron spectra in 6 = collisions and all 
spectra for 6^0 collisions (including all flow anisotropies such as vi which vanish at 6 = 0) are 
then parameter-free predictions of the model. Note that all calculated hadron spectra include 
feeddown from decays of unstable hadron resonances. 

4. Successes of ideal fluid dynamics at RHIC 

4-1- Hadron momentum spectra and radial flow 

Figure El shows the single particle pj_-spectra for pions, kaons and antiprotons (left panel) 
as well as O baryons (right panel) measured in Au+Au collisions at RHIC, together with 
hydro dynamical results ^2j. In order to illustrate the effect of additional radial flow generated in 
the late hadronic stage below T cr , two sets of curves are shown: the lower (blue) bands correspond 
to kinetic decoupling at T cr = 165 MeV, whereas the upper (red) bands assume decoupling at 
Tdec = 100 MeV. The width of the bands indicates the sensitivity of the calculated spectra to 




Figure 2. Negative pion, kaon, antiproton, and spectra from central Au+Au collisions at 
V / s = 200AGeV measured at RHIC 16 . The curves show hydrodynamical calculations ^2] (see text). 



an initial transverse flow of the fireball already at the time of thermalization ^2]. In the 
hydro dynamic simulation it takes about 9-10 fm/c until most of the fireball becomes sufficiently 
dilute to convert to hadronic matter and another 7-8 fm/c to fully decouple [£]. Figure |2] shows 
that by the time of hadronization the dynamics has not yet generated enough radial flow to 
reproduce the measured p and f2 spectra; these heavy hadrons, which are particularly sensitive to 
radial flow, require the additional collective "push" created by resonant quasi-elastic interactions 
during the fairly long-lived hadronic rescattering stage. The flattening of the p spectra by radial 
flow provides a natural explanation for the (initially puzzling) experimental observation that for 
p± > 2 GeV/c antiprotons become more abundant than pions |15j . 

As shown elsewhere (see Fig. 1 in |17j ) , the model describes these and all other hadron spectra 
not only in central, but also in peripheral collisions, up to impact parameters of about 10 fm, 
and with similar quality. No additional parameters enter at non-zero impact parameter. 

4-2. Elliptic flow 

Figure 01 shows the predictions for the elliptic flow coefficient V2 from Au+Au collisions at RHIC, 
together with the data |18l 119] . For impact parameters b < 7fm (n c h/n max > 0.5) and transverse 
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Figure 3. Left: p^-averaged elliptic flow for all charged hadrons from 130 AGeV Au+Au collisions, as 
a function of collision centrality (n c h is the charged multiplicity at y = 0). The curves are hydrodynamic 
calculations with different choices for the initial energy density profile (see [23j). Right: Differential 
elliptic flow V2(px) for identified hadrons from minimum bias Au+Au collisions at 200 A GeV |18U19ll2*T] . 
together with hydrodynamic curves from |22| . 

momenta p± < 1.5GeV/c the data are seen (left panel of Fig. OJ) to exhaust the upper limit 
for V2 obtained from the hydrodynamic calculations. For larger impact parameters b> 7fm the 
p_i_-averaged elliptic flow V2 increasingly lags behind the hydrodynamic prediction; this will be 
discussed in detail in Sec. 15.21 As a function of p± (right panel of Fig. I2J) the elliptic flow of all 
hadrons measured so far is very well described by hydrodynamics, for p± < 1-5 — 2GeV/c. In 
particular the hydrodynamically predicted mass splitting of v<i at low p± is perfectly reproduced 
by the data. This mass splitting depends on the EOS [Hj, and the EOS including a quark- 
hadron phase transition used here describes the data better than one without phase transition 
(see Fig. 2 in J7]). Ideal fluid dynamics with a QGP EOS thus gives an excellent and very 
detailed description of all hadron spectra below pj_ = 1.5 GeV/c. Since this pj_-range includes 
more than 99% of all produced hadrons, it is fair to say that the bulk of the fireball matter formed 
in Au+Au collisions at RHIC behaves very much like a perfect fluid. 

4-3. Final source eccentricity in coordinate space 

While spectra and elliptic flow reflect the momentum structure of the hadron emitting source, 
its spatial deformation can be tested with 2-pion Hanbury Brown- Twiss (HBT) correlations 
measured as a function of the azimuthal emission angle |231 124j . Even though the initial spatial 



deformation of the reaction zone in non-central Au+Au collisions at RHIC is finally completely 
gone, many pions are already emitted from earlier times when the spatial deformation is still 
significant. For Au+Au collisions at 6 = 7fm, the spatial eccentricity of the time-integrated 
hydrodynamic pion emission function is (e x ) = 0.14 (or 56% of its initial value £ x {to) = 0.25) |23j . 
Using azimuthally sensitive pion HBT measurements, the STAR Collaboration has measured in 
the corresponding impact parameter bin [2H] (e x ) = 0.11 ± 0.035 (or 45 ± 15% of the initial 
deformation). This agreement can be counted as another success for the hydrodynamic model. 

5. Viscous effects at RHIC 

5.1. QGP viscosity 

As evident in the right panel of Fig. [31 the hydrodynamic prediction for V2(p±) gradually 
breaks down above p± > 1.5GeV/c for mesons and above p± > 2.2GeV/c for baryons. The 
empirical fact [2J that both the pj_-values, where this break from hydrodynamics sets in, and 
the saturation values for V2 at high p± for baryons and mesons are always (i.e. for all collision 
centralities!) related by 3 : 2 (i.e. by their ratios of valence quark numbers), independent of 
their masses, tells its own interesting story (see e.g. [213 ) : It strongly suggests that in this 
p± region hadrons are formed by coalescence of color-deconfined quarks, and that the elliptic 
flow is of partonic origin (i.e. generated before hadronization), with a p^-shape that follows 
hydrodynamics at low p± up to about 750 MeV and then gradually breaks away |2(ij . 

Since i>2(p±) is a measure for the relatively small differences between the in-plane and out- 
of-plane slopes of the p± spectra, it is more sensitive to deviations from ideal fluid dynamic 
behaviour than the angle-averaged slopes. Two model studies (23 I2B] showed that V2 reacts 
particularly strongly to shear viscosity. As the mean free path of the plasma constituents (and 
thus the fluid's viscosity) goes to zero, V2 approaches the ideal fluid limit from below [22] (see 
Fig-Ek) • At higher transverse momenta it does so more slowly than at low p± [23] , approaching 
a constant saturation value at high p±. The increasing deviation from the ideal fluid limit for 
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Figure 4. Left: Elliptic flow from a parton cascade [2Hj, compared with STAR data, for different 
parton-parton scattering cross sections. Larger cross sections lead to smaller mean free paths. Right: 
Pcrturbative effects of shear viscosity on the elliptic flow V2(p±) [2Hj (see text). 

growing p^_ is qualitatively consistent with the expected influence of shear viscosity: Teaney |28| 
showed that the lowest order viscous corrections to the local equilibrium distribution increase 
quadratically with p± so that v-i remains increasingly below the ideal fluid limit as p± grows 
(see Fig. Eh). From Fig. 2b Teaney concluded that at RHIC the normalized sound attenuation 
length ^f- = 37yf (where r\ is the shear viscosity, T the temperature and s the entropy density) 
cannot be much larger than about 0.1. This puts a stringent limit on the dimensionless ratio rj/s, 
bringing it close to the recently conjectured absolute lower limit for the viscosity of rj/s = fij (47r) 
|3()j . This would make the quark-gluon plasma the most ideal fluid ever observed |3()j . 

These arguments show that deviations from ideal fluid dynamics at high p± must be expected, 
and that they can be large even for fluids with very low viscosity. At which p± non-ideal 



effects begin to become visible in V2(p±_) can be taken as a measure for the fluid's viscosity. To 
answer the quantitative question what the RHIC data on partonic elliptic flow and its increasing 
deviation from ideal fluid behaviour above p± > 750MeV/c imply for the value of the QGP 
shear viscosity r] requires a numerical algorithm for solving viscous relativistic hydrodynamics 
(see jSJ f° r an update). 

5. 2. Viscosity of the hadron resonance gas 

Ideal fluid dynamics also fails to describe the elliptic flow V2 in more peripheral Au+Au collisions 
at RHIC and in central and peripheral collisions at lower energies (see Fig. Et), as well as at 
forward rapidities in minimum bias Au+Au collisions at RHIC |H2] . Whereas hydrodynamics 
predicts a non-monotonic beam energy dependence of V2 (Fig. e 6„), with largest values at 
upper AGS and lower SPS energies, somewhat lower values at RHIC and again larger values at 
the LHC, the data seem to increase monotonically with y/s. 
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Figure 5. Left: Scaled elliptic flow v^/e (where e denotes the initial spatial eccentricity) vs. the charged 
multiplicity density per unit initial transverse overlap area S Right: Elliptic flow V2 for minimum 
bias Au+Au collisions at various collision energies, parametrized by the final charged multiplicity density 
at midrapidity See text for discussion. 

Such a monotonic rise is consistent with "hybrid" calculations by Teaney |14j (Fig. JSJa) where 
the fireball undergoes ideal fluid dynamic evolution only while in the QGP, followed by hadronic 
kinetic evolution using RQMD after hadronization. Figured shows several curves corresponding 
to different equations of state during the hydrodynamic evolution (see ^U), with LH8 being 
closest to the data. The difference between the points labelled LH8 and the hydrodynamic 
values at the top of the figure is due to the different evolution during the late hadronic stage. 
Obviously, at lower collision energies and for impact parameters b ~ 7 fm (which dominate V2 in 
minimum bias collisions), ideal fluid dynamics continues to build additional elliptic flow during 
the hadronic stage, but RQMD does not. The initial energy densities are smaller than at RHIC 
and the fireball does not spend enough time in the QGP phase for the spatial eccentricity e to 
fully disappear before entering the hadron resonance gas phase. Anisotropic pressure gradients 
thus still exist in the hadronic phase, and ideal fluid dynamics reacts to them according to the 
stiffness of the hadron resonance gas EOS (p ~ 0.15e). Teaney 's calculations |14j show that 
RQMD responds to these remaining anisotropies much more weakly, building very little if any 
additional elliptic flow during the hadronic phase. The hadron resonance gas is a highly viscous 
medium, unable to maintain local thermal equilibrium. The failure of the hydrodynamic model 
in situations where the initial energy density is less than about 10GeV/fm 3 [21] is therefore 
likely not the result of viscous effects in the early QGP fluid, but rather caused by the highly 
viscous late hadronic stage which is unable to efficiently respond to any remaining spatial fireball 
eccentricity. Similar arguments hold at forward rapidities at RHIC |H4| where the initial energy 
densities are also significantly smaller than at midrapidity while the initial spatial eccentricities 
are similar. 



The large hadron gas viscosity spoils one of the clearest experimental signatures for the 
quark-hadron phase transition, the predicted |Hj non-monotonic beam energy dependence of V2 
which was already described in Sec. [2] and is shown in Fig. EH- As one comes down from infinite 
beam energy, vi is predicted to first decrease (due to the softening of the EOS in the phase 
transition region) and then recover somewhat in the moderately stiff hadron gas phase. The 
hadron gas viscosity spoils this recovery, leading to an apparently monotonous decrease of vi 
with falling beam energy (Fig. [5J). However, recent PHENIX data from Au+Au collisions at 
yfs = 62 A GeV (33] indicate that this decrease may not be quite as monotonous as suggested by 
Fig. Figure shows that, at fixed p±, the elliptic flow V2(p±) is essentially constant over 
the entire energy range explored at RHIC (from 62 to 200 A GeV), decreasing only when going 
further down to SPS and AGS energies. When integrated over p±, this turns into a monotonous 
behaviour as in Fig. due to the steepening p± spectra at lower y/s. While Fig. 03 does not 
confirm the hydro dynamically predicted rise of v% towards lower y/s, it may still reflect this 
predicted non-monotonic structure in the elliptic flow excitation function, after strong dilution 
by hadronic viscosity effects. Obviously, many and more systematic hybrid calculations of the 
type pioneered by Teaney |14j are necessary to explore to what extent we can eventually prove 
the existence of a QCD phase transition using elliptic flow data in the SPS-RHIC energy domain. 
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Figure 6. Left: Excitation function of radial and pj_-integrated elliptic flow for b = 7fm Pb+Pb or 
Au+Au collisions The horizontal axis gives dN 7r /dy(b=7im), and the horizontal arrows indicate how 
dNir/dy was expected to correspond to the beam energy ranges covered by SPS, RHIC and LHC before 
the RHIC data were available. Right: Elliptic flow at fixed p± = 0.65 and 1.75 GeV/c from A+A collisions 
with A w 200 at a variety of beam energies [HE] . 

Obviously, it will be important to confirm non-viscous fluid behaviour at higher initial energy 
densities than so far explored, by extending Fig. |HK to the right and verifying that vi/e settles 
on the hydrodynamic curve. This can be done with Pb+Pb collisions at the LHC, or with full- 
overlap U+U collisions at RHIC [I]. In addition, the large spatial source deformations achievable 
in full-overlap side-on-side U+U collisions allow for decisive systematic studies of the non-linear 
path-length dependence of QCD radiative energy loss of fast partons. In [3] we give quantitative 
arguments why a U+U collision program should be seriously considered at RHIC. 

6. Conclusions 

The collective flow patterns observed at RHIC provide strong evidence for fast thermalization 
at less than 1 fm/c after impact and at energy densities more than an order of magnitude above 
the critical value for color deconfinement. The thus created thermalized QGP is a strongly 
coupled plasma which behaves like an almost ideal fluid. These features are first brought out in 
heavy-ion collisions at RHIC because only there the initial energy densities and QGP life times 
are large enough for the ideal fluid character of the QGP to really manifest itself, in the form 
of fully saturated hydrodynamic elliptic flow, undiluted by late non-equilibrium effects from the 
highly viscous hadron resonance gas which dominates the expansion at lower energies. 



We are now ready for a systematic experimental and theoretical program to quantitatively 
extract the EOS, thermalization time and transport properties of QGP and hot hadronic matter. 
This requires more statistics and a wider systematic range for soft hadron production data, but 
more importantly a wide range of systematic simulation studies with the "hydro-hadro" hybrid 
algorithms and, above all, a (3+l)-dimensional viscous relativistic hydrodynamic code (see [HJ 
for more on that). 
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